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Abstract 



The infrared behaviour of vertex functions in an SU{N) Yang-Mills theory in Lan- 
dau gauge is investigated employing a skeleton expansion of the Dyson-Schwinger 
equations. The three- and four-gluon vertices become singular if and only if all ex- 
CN ■ ternal momenta vanish while the dressing of the ghost-gluon vertex remains finite 

I in this limit. The running coupling as extracted from either of these vertex func- 

' tions possesses an infrared fixed point. In general, diagrams including ghost-loops 

. dominate in the infrared over purely gluonic ones. 

(N 
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Fifty years after the formulation of Yang-Mills theory its infrared (IR) struc- 
ture is still largely unknown despite the fact that this knowledge is central to 



;h ' any effort in understanding the strong interactions from first principles. It has 

^ long been conjectured that IR enhancements are present. Indeed they are nec- 

essary to explain confinement. Such IR enhancements may also be the reason 
that no explicit glue becomes visible in the low-lying hadron mass spectrum. 

Lattice calculations include in principle all non-perturbative features of Yang- 
Mills theories but are in practice limited by the finite lattice volume in the 
study of possible IR singularities [1,2,3,4]. A complementary non-perturbative 
continuum method is provided by the Dyson-Schwinger equations (DSEs). 
In Landau gauge the DSEs for the ghost and gluon propagators have been 
analytically solved in the IR assuming ghost dominance [5,6,7,8,9,10,11]. Lat- 
tice and DSEs are complementary and yet they agree on the propagators' IR 
behaviour: there is clear evidence for an IR finite or even vanishing gluon prop- 
agator and a strongly diverging ghost propagator, in accordance with both, 
the Kugo-Ojima confinement criterion [12] and the Gribov-Zwanziger scenario 
[13]. 
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The running of the gauge couphng is intimately related to the momentum 
dependence of the primitively divergent vertex functions in an SU (A^c) Yang- 
Mills theory. For sufficiently large momenta the coupling can be calculated 
from perturbative corrections to these vertex functions [14]. The breakdown 
of perturbation theory is signaled by Landau poles, unphysical singularities at 
non- vanishing space-like momenta. By simply imposing analyticity for space- 
like momenta extrapolations to the IR have been performed in so-called An- 
alytic Perturbation Theory [15,16,17]. Typically these studies find a well- 
behaved coupling at all momenta and an IR fixed point. Qualitatively this 
agrees with the findings from the genuinely non-perturbative continuum ap- 
proaches. In the latter the IR behavior of the ghost-gluon interaction in Lan- 
dau gauge has been determined either from DSEs or the Exact Renormaliza- 
tion Group Equations (ERGEs) [8,9,18,19] and yield an IR fixed point with 
a{0) ~ 8.9/iVc. The corresponding couplings from the three- and four-gluon 
vertex functions have not yet been studied with these techniques. Within per- 
turbation theory the principle of gauge invariance leads to the universality of 
the gauge coupling: All three couplings are equal up to the well-known renor- 
malization scheme dependencies [14,20,21] at every finite order. Although a 
general proof is lacking one expects such relations to also hold in the full 
non-perturbative theory. 

In this letter we propose a method to investigate the IR behaviour of Greens 
functions with an arbitrary number of ghosts and gluons. We detail the results 
for the gluon-self interaction vertices thus completing our knowledge about the 
primitively divergent vertex functions {i.e. those appearing in the renormal- 
ized Lagrangian) and the corresponding running coupling(s) in the deep IR. 
We construct a skeleton expansion for each DSE employing the fully dressed 
primitively divergent n-point functions of Yang- Mills theory. To all orders in 
this expansion the three- and four-gluon vertex functions are IR singular if 
and only if all external momenta vanish. The corresponding IR exponents are 
hereby proportional to the IR exponent of the ghost renormalization function 
such that the corresponding couplings have IR fixed points. The results pre- 
sented in this letter verify the self-consistency of ghost dominance in a quite 
general way and prove that ghost dominance is a successful guiding princi- 
ple when determining the IR behaviour of Yang-Mills Greens functions in the 
confining phase. 

In Landau gauge, the ghost and gluon propagators in Euclidean momentum 
space are described by the renormahzation functions G(p^) and Z{p'^): 



On the other hand, the three- and four-point functions feature tensor struc- 
tures not present in the bare Lagrangian thus requiring multiple such scalar 
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functions. We will not specify these yet and discuss first the behaviour of the 
scalar amplitude multiplying only the tensor defined on the level of the bare 
Lagrangian. E.g. for the three-gluon vertex this tensor reads 

<^WM2(Pl - P2)m3 + <^M2M3(P2 - P3)/.i + (^M3w(P3 " Pi) 1x2 ■ (2) 

In a setting where all external momenta (pi)^ vanish as a single momentum 
scale — > this tensor will be multiplied by an amplitude Hl^{p^) in the fully 
dressed three-gluon vertex. For the four-gluon and ghost-gluon vertex func- 
tions the amplitudes Hi^{p^) and Hf^~^{p^) are defined analogously. Starting 
from the established result [8,9,18,19] 

Z(p')^(pY\ G(p')^(pY\ (3) 

we will show in the following that for — > 

H',%p') - {pY'^ , Ht%p') {pY'^ , Ht-%p') const. (4) 

where the parameter k is determined from the propagator DSEs, see below. 
Eqs. (3) imply that the minimum of a gluon's dispersion relation does not occur 
at zero momentum. ^ In addition, it concurs with Zwanziger's picture [13]: the 
geometric degrees of freedom dominate IR Yang-Mills theory. Eqs. (4), on the 
other hand, state that the self-interactions of gluons become large at very low 
momentum. The connected three- and four-gluon functions hereby violate the 
bound conjectured in ref. [23], although only at one specific kinematical point. 

It will be demonstrated in the following that each of the eqs. (3) and (4) is 
correct when the others arc assumed as given. At the end, one thus obtains 
a self-consistent solution of the DSE system in the IR. The starting point is 
provided by the observation that in Landau gauge the ghost-gluon vertex is 
UV finite and remains bare for vanishing incoming ghost momentum [24,25]. 
This is easily shown using the DSE for the full ghost-gluon vertex, see fig. 1. 
The bare ghost-gluon vertex in the interaction diagram is proportional to the 
internal loop momenum l^. Due to the transversality of the gluon propagator 
Df^,^ one has l^D^y{l—q) = q^D^yil—q). Thus, the interaction diagram vanishes 
for — > 0. This argument would only be invalidated if the two-ghost-two- 
gluon scattering kernel had an IR divergence. However, recent lattice results 
and numerical DSE studies of the ghost-gluon vertex [26,27] agree with a 
bare vertex in the IR. From our analytical study it also follows that such a 
divergence is absent. 

The IR behaviour of two-point functions has been adressed within ERGEs 
and DSEs [6,7,8,9,10,18,19]. We briefiy discuss the IR aspects of the ghost 
propagator DSE depicted in fig. 2. With an IR finite ghost-gluon vertex, see 

^ The phenomenological consequences of an IR suppressed gluon propagator can 
for example be infered from ref. [22]. 
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Fig. 1. Ghost-gluon vertex DSE. 




Fig. 2. Ghost propagator DSE. 



the last of eqs. (4), and a power law ansatz for the dressing functions at low 
external p^, Z{p'^) = A (p^)" and G{p'^) = B (p^)^, a self-consistent solution 
arises with k := —(3 = 2a from matching the exponents of both sides of 
the ghost DSE [28]. Hereby the formula for a scalar integral in d dimensions, 
depending on one external momentum 



J dq{q} {{q-p) ) - [p ) n^T{-a)T{-b)Tid + a + b) 

has been used. The actual value of k follows from IR consistency of the ghost 
with the gluon propagator DSE, and depends slightly on the details of the 
ghost-gluon vertex' (finite) dressing [9]. For a bare vertex one obtains k, = 
(93 — vl201)/98 ~ 0.595. Note, however, that the precise value of k, will be 
irrelevant for all arguments presented in the following as long as < k < 1 
[9,28] and k^1/2[9]. 

In a first step we establish the three-gluon vertex power law in eq. (4). The 
corresponding DSE (see fig. 3) includes four- and five-point functions which 
we treat in a skeleton expansion in terms of the fully dressed, primitively 
divergent Greens functions. Furthermore, we develop counting rules for the 
IR degree of divergence of a diagram. To lowest order in the expansion we can 
identify each subdiagram with its leading IR singularity. We will demonstrate 
that all higher order diagrams are either equally or less IR singular. If the three 
momenta entering the vertex are much smaller than Aqcd the integral will be 
dominated by loop momenta also smaller than Aqcd due to the denominators 
of the propagators. The dressing functions of these propagators are then the 
simple power laws given in eqs. (3), allowing to integrate the loops analytically. 
The ghost-loop diagram (a) in fig. 3 is IR leading. We thus evaluate this 
diagram first and then substitute the resulting power law for the three-gluon 
vertex into the other diagrams to check for self-consistency. This is simplest 
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Fig. 3. Exact Dyson-Schwinger equation for the three-gluon vertex and lowest order 
in a skeleton expansion of the four- and five-point functions. All internal propagators 
in the diagrams are to be understood as fully dressed. 

demonstrated at the symmetric kinematical point 

pI=pI=pI = P^ (Pi • P2) = (pi • Ps) = {P2 ■ Pa) = -^P^- (6) 

The tensor structure of the three-gluon vertex then reduces from fourteen 
tensors at general kinematics to three independent ones [14]. Omitting the 
overall color factor and the delta function expressing momentum conservation 
the vertex is then given by 



(Pi , P2, Ps) = hI^{p^) {5 1X11X2 (Pi - P2)m3 + cychc permutations) 
-HI\p^)Ip^ (P2 -P3)mi(P3 -Pi)m2(Pi -P2)m3 
+ hI'{p')/p' {pi,^P2,,Pa,, - Pi,,P2,,Pz,,) (7) 

with the momentum routing and indices as given in fig. 4. 

We are interested in the limit — > where the ghost-gluon vertex becomes 
bare again [24], T^^{p,q) = iq^ with being the momentum of the outgoing 
ghost. Substituting the power laws, eqs. (3), for the ghost and gluon dressing 
functions and employing (5) the IR behavior of the dressing functions if^^, 
and i?3^ of the vertex is determined analytically to be 
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Fig. 4. Momentum routing in the ghost-loop diagram of the three-gluon vertex. All 
external momenta flow into the loop. Internal propagators are fully dressed. 
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(9) 
(10) 

lO 

with 

^ ' 327r2 r(l + fi;)r(2 + 2K)r(3-3K) ^ ' 

and B being the leading IR coefficient of the ghost. Crucial to these expressions 
are the momentum power laws. All three dressing functions are proportional to 
j^p2^(-3K)^ z.e. with K > [28] the three-gluon vertex is IR singular, and the first 
of eqs. (4) is established. Note that the degree of singularity is the same for all 
three dressing functions, and given by three times — k, the latter being the IR 
exponent of the ghost propagator. In the following it will become evident that 
the IR anomalous dimension of a diagram is simply the sum of IR exponents 
of propagators and vertices constituting the diagram. All trivial dimensions 
of the vertex are due to the tensor structure whereas the dressing functions 
iff^, and H^^ are dimensionless. Any anomalous dimension appearing in 
the diagram of fig. 4 via nontrivial vertex or propagator dressings therefore 
has to appear in the functions H^^, ifg^ and -^3^. Therefore we can determine 
the degree of divergence of Greens functions just by counting IR anomalous 
dimensions. As is the only scale, {p^Y with p = —3k has necessarily to arise 
in this case. 



Note that a dressed ghost-gluon vertex would not have changed the running of 
H^^, and H^^ with momentum, since in Landau gauge it carries no overall 
anomalous dimension. IR corrections to the vertex dressings can therefore be 
typically expressed as quotients [9]. Such dressings cannot change the 

momentum dependence in eqs. (8-10). The following observation is important: 
The diagram in fig. 4 becomes IR singular if and only if all three incoming 
momenta (pi)^, {P2Y (Ps)^ approach zero. This can be proven in general 
kinematics.^ A convenient kinematic section to see this directly is e.g. pro- 

2 Using appropriate projectors to single out the scalar dressing functions of the 
vertex one always ends up with scalar 'massless' triangle integrals with nontrivial 
powers of internal momenta. These integrals can be performed along the methods 
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Fig. 5. Typical terms of the skeleton expansion of the ghost-gluon scattering kernel 
in the DSE for the three-gluon vertex. Internal propagators are fully dressed. 

vided by pi — constant, p2 — bpi with b a constant that we will allow to 
approach zero. Due to momentum conservation the three momenta are then 
collinear and only one of them is small. With a convenient choice for the tensor 
basis, the leading coefficient calculated from the ghost triangle graph is 



In the limit 6 — > the integrand is singular but integrable, and thus the 
corresponding diagram yields a finite contribution. 

In the DSE for the three-gluon vertex no other term in the skeleton expan- 
sion has a larger power IR divergence. To show this we count powers for the 
remaining diagrams in fig. 3 and compare the result with eqs. (8-10). Employ- 
ing eqs. (3-4) and denoting the remaining four diagrams as (6) — (e), c.f. fig. 
3, we obtain their IR anomalous dimensions by summing up the anomalous 
dimensions of all propagators and vertex functions of a diagram. Hereby a 
gluon propagator contributes an anomalous dimension of 2k. a ghost propa- 
gator one of (— k), the three-gluon vertex one of (— 3/t) and the ghost-gluon 
vertex carrys zero anomalous dimensions, c.f. eqs. (3) and (4). This yields for 
their respective IR exponents p: 

(b) p = 3 X 2k, + 2 X (-3k) = 0, (c) p = 5 x 2k + 3 x (-3k) = k, 

(d) p = 2 x 2k + 1 x (-3k) = k, (e) p = 2 X 2k + 1 X (-3k) = k. (13) 

Thus they are all subleading if compared to the ghost triangle diagram (a). 

Finally we investigate further terms from the skeleton expansion. Typical cor- 
rections to the ghost-triangle diagram are given in fig. 5. The only scale ap- 
pearing in the symmetric momentum kinematics is (p^), thus we can determine 
the degree of divergence of the diagrams by counting anomalous dimensions. 




/ 



d^k b{b - 3) {kYi -k-pl) + {b- l)2ek ■ pi 
(27r)4 (62 -36 + 3) 



G{k') G{{k + bp,f) G{{k - p,)') 
k^ {k + bpif {k-pif 



(12) 



described in refs. [30,31]. As we have checked explicitly, no IR singularities arise 
unless all three squared momenta vanish. 
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Fig. 6. Insertions generating the loop expansion in fig. 5. One should also allow for 
the possibility of upgrading a three to a four gluon vertex. 

Denoting the higher order diagrams by (b) — (d) we obtain for their IR anoma- 
lous dimensions 



with no contribution from the ghost-gluon vertex, as discussed above. Thus all 
these diagrams contribute to the same IR order as the ghost-loop diagram (a). 

Arbitrary diagrams in the expansion can be generated with the four elements 
in fig. 6, by repeated insertion into a given lowest-order diagram. According to 
our counting rules based on eqs. (3) and (4) these insertions carry an overall 
zero IR anomalous dimension. Therefore in general, starting from a given dia- 
gram in the skeleton expansion with an IR anomalous dimension p all possible 
insertions generate diagrams of the next order in the skeleton expansion with 
the same degree of divergence. By induction this is true for the whole series 
of skeleton diagrams up to any order. Therefore the degree of IR divergence of 
a fully dressed Greens function in the presence of only one external scale can 
be already read off by counting the anomalous dimensions of the lowest order 
diagrams in the skeleton expansion. This general rule especially implies that 
higher order corrections to the diagrams (6) — (e) of fig. 3 are all IR subleading. 
Thus the IR exponent p = —3k for the three-gluon vertex is established. ^ 



Fig. 7. First order expansion of the DSE for the four-gluon correlation. Internal 
propagators are to be understood as fully dressed. 

The lowest order terms in the skeleton expansion of the DSE for the four-gluon 
vertex are given in fig. 7. Here we consider the special momentum configuration 
where all four momenta flowing into the loop are equal in magnitude and 

^ This result has been anticipated in [32] employing the STI Zi{iJ.) = Z^{p) / Z^{ix). 
In the presence of only one momentum scale all renormalized dressing functions run 
with where a is the corresponding anomalous dimension. Therefore the 

IR powers in Eq.(3) imply that Zi{p) ~ (fp)^'^ and p = —3k follows directly. 



(6) p= (3-2 + 3-(-l) + 2-(-3))fi: 
(c) p = (2-2 + 7- = -3k, 

{d) p= (4-2 + 8-(-1) + 1-(-3))k 



—3k, 



—3k, 



(14) 




pairwisc parallel or antiparallcl. Denoting the ingoing momenta with Pi=i,,A we 
have Pi=P2=P3=p'i^ P^ and all three Mandelstam variables s = (^1+^2)^; 
t = ip2 H-Ps)^ and u = (pi + Ps)^ are either zero or proportional to p^. With 
the IR counting rules established above we can then determine the degree of 
divergence of the three one- loop diagrams (a), (b) and (c). We obtain 

(a) p = -4k, (6) p = (4 • 2 + 3 • (-3)) k = -k, 

(c) p = (3 • 2 + 1 • (-3) + 1 • (-4)) K = -K. (15) 

Here the ghost box diagram is IR dominating, and yields —Ak as IR exponent 
for the four-gluon vertex. 

Finally, we have to check whether eqs. (3) remain valid assuming the new 
results for the gluon-selfinteraction vertices. We therefore reconsider the DSE 
for the gluon propagator given diagramatically in fig. 8. Again due to the 
denominators of the propagators all diagrams are dominated by loop momenta 
similar to the external scale p^. We can thus determine the IR behaviour for 
each diagram by the counting rules developed above. Note that the tadpole 
diagram contains only the bare four-gluon vertex and is independent of the 
external momentum. ^ For the other diagrams we obtain the following degree 
of IR divergence: 

(a) p = -2k, (b) p = K, (c) p = 2k, (d) p = 2k. (16) 

Clearly, the ghost-loop diagram is the leading diagram in the IR, and the 
analysis leading to eqs. (3) is thus not altered. 

With analogous techniques one can straightforwardly prove that the generic IR 
behaviour of a 2n-ghost-m-gluon amputated and connected Greens function 
is given by (p2-)(n-m)K; 

To summarize this part, assuming that skeleton expansions of higher n-point 
Greens functions are well-defined we have shown: 

• The degree of nonperturbative IR divergence of an n-point function is given 
already by the highest degree of divergence present in the lowest order skeleton 

expansion of the diagrams appearing in its DSE. 

• In Landau gauge Yang-Mills theory the three-gluon vertex behaves as (p^)"^" 
and the four-gluon vertex as (p^)"^** in the IR. 

A consistency check on this IR behaviour is, of course, given by its implication 
on the IR value of the running coupling as infered from these vertex functions. 
As we will see, these results allow for universality of the IR fixed point in Yang- 
Mills theory. A nonperturbative expression for the running coupling has been 



^ The tadpole diagram can even be completely absorbed in the renormalization. 
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Fig. 8. The Dyson- Schwinger equation for the gluon propagator. 

derived in the context of DSEs from the ghost-gluon vertex [6]. The Slavnov- 
Taylor identity (STI) 

Z^ = Z,Z^Zy\ (17) 

relates the vertex renormahzation factor Z^ with the corresponding factor 
Zg for the couphng g and the ghost and gluon fields, Z^^ ^ Z^^^. With the 
definition a — g'^/{A7r) and the relation g'(A^) = A^)g'(//^) one obtains 

for the renormalized coupling at the renormahzation point /i^. The bare cou- 
pling a(A^) depends on a regulator A^. Since the ghost-gluon vertex is finite 
in Landau gauge, one may choose Z\ = \. The bare and renormalized ghost 
and gluon dressing functions are related by 

Go(p^ A^) = G{v\ //')^3(//^ A^) , zo(p^ a^) = z{v\ /^')^3(/^^ a^) . (19) 

Substituting these relations into eq.(18) we obtain 

a{ix^)G\v\ i^^)Z{v\ i^^) = a(A^)G^(p', A^)^o(p^ A^) ■ (20) 

Note that the r.h.s. is independent of [j? and thus the l.h.s. is renormahzation 
group invariant [29]. Evaluating the l.h.s. once at an arbitrary renormahzation 
point //^ and once at — we obtain 

a{p■')^a{^J?)G\p^^^'')Z{p\^^'') (21) 

where we have exploited the renormahzation condition G"^ {p^ , p^) Z {p"^ , p'^) = 1 
for the renormalized dressing functions. Evaluating Q!(p^) in the ultraviolet one 
recovers the well known perturbative coupling in the MOM-scheme which can 
be related to the MS'-couphng by standard techniques [14] . 

The IR behaviour of this coupling can be read off from the power law behaviour 
of the ghost and gluon dressing functions in eqs. (3). Due to the interrelated 
exponents in these expressions the IR momentum dependence cancels and 
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leads to a fixed point at = 0. The precise value of a{0) depends on k. 
It has been calculated in both, DSEs and ERGEs, employing several kinds 
of truncation schemes [6,7,8,9,10,18,19]. Neglecting IR corrections from the 
dressing of the ghost-gluon vertex yields 



For the running coupling based on the three-gluon vertex we employ the STI 

Z, = Zg ZT (23) 

relating the renormalization factor of the three-gluon vertex to the one of the 
couphng and the gluon fields. The vertex renormalization factor Z\ relates the 
regularized and the renormalized three-gluon vertex via 

r'^^'(pi,P2,P3,/^') = rr(pi,P2,P3, A') A') . (24) 

Evaluating the vertex at the symmetric point Pi = pi — Pa = obtains 

El\v\ 1^') = Hl\p\ A') Zl(/.^ A^) (25) 

and thus a relation between the renormalized and the regularized H^^ as de- 
fined in eq. (7). One can also derive an expression for the three-gluon coupling: 

= «o'(^') {H!r{p',f^') Z\p\^i') (26) 
where we exploited the renormalization condition (iff^)^(p^,p^)Z^(p^,p^) = 1. 
This condition involves only H^^ because it multiplies the tensor structure of 
the three-gluon vertex containing the primitive UV divergence. Again the r.h.s. 
of this equation denotes an RG invariant quantity. This coupling is related to 
the one of the ghost-gluon vertex by a known scale transformation [14]. 

RecaUing that Hl^{p^) ~ {p^Y^'' and Z(p2) ^ (p2)2« ohi^m 

a^%p^ ^ 0) - const. /Nc , (27) 

i.e. the running coupling taken from the three-gluon vertex possesses an IR 
fixed point in accordance with the coupling determined from the ghost-gluon 
vertex. The explicit l/A^^c-dependence stems from the fact that ~ l/A'c 
and the functions and Z are independendent of N^, as can be seen from 
their DSEs. Note that also other renormalization conditions for the three- 
gluon vertex involving and H^^ are possible [14]. The expression (26) for 
the coupling would then include a corresponding linear combination of the 
functions H^^ , and H^^ . All of these functions are proportional to (p^)"^" 
in the IR, therefore such a redefinition of the coupling would not affect the 
presence of the IR fixed point. 



11 



Finally we examine the four-gluon coupling. Its STI is — Z'^Z\^ and thus 
we obtain the relation 

E{\V\ /^') = ii\\v\ A') ^4(/^^ A') . (28) 
The running coupling from the four-gluon vertex is therefore given by 

a^^(p^) = a'^{ii^) E\\v\ /^') /.^) (29) 

where analogously to the three-gluon vertex case the renormalization condition 
H^^ {p'^ , p'^) Z'^ {p'^ , p^) = 1 has been employed. Recalling that i/f^(p^) ~ (p^)"^** 
and Z^p"^) ~ (p^)^'^ in the IR we again obtain an IR fixed point, 

a'^^p^ ^ 0) ~ const. /Nc. (30) 



To summarize, we have demonstrated on a qualitative level that the running 
couplings from the ghost-gluon, three-gluon and four-gluon vertices have a 
universal, nontrivial fixed point in the IR. The value of the fixed point for 
the ghost-gluon coupling is known to depend slightly on the IR dressing of 
the ghost-gluon vertex and falls into a small window 2.5 < a{0) < 3 for 
A*",, = 3 [9]. The corresponding value for the other couplings can be calculated 
from the vertex-DSEs. In the skeleton expansion the dressed three- and four- 
gluon vertices receive their IR leading contributions from an infinite number 
of diagrams belonging to a certain subclass. These can be constructed from 
the ghost-triangle diagram by repeated insertion of the diagrammatic pieces 
given in fig. 6. In the IR, all these diagrams have the same dependence on 
one external momentum scale. In principle, the momentum independent coef- 
ficients can be calculated order by order in the skeleton expansion, although 
its convergence properties are yet to be determined. The nonperturbative ex- 
pressions for the three- and four-gluon vertices are IR singular provided all 
external squared momenta approach zero. This singularity, however, is not 
strong enough to compensate the zeroes in the gluon dressing functions of 
attached gluon propagators. Thus in the skeleton expansion of all DSEs the 
ghost-loop diagrams are leading in the IR. This is in agreement with a picture 
recently advocated by Zwanziger [13]: the geometric degrees of freedom, i.e. 
the Faddeev-Popov determinant, dominate IR Yang-Mills theory. 
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